Cosmology with space-based gravitational-wave detectors 
— dark energy and primordial gravitational waves — 
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Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect 
^ 10^ neutron-star (NS) binaries and determine the luminosity distances to the binaries with high 
precision. Combining the luminosity distances with cosmologically-induced phase corrections on 
the GWs, cosmological expansion out to high redshift can be measured without the redshift deter- 
minations of host galaxies by electromagnetic observation and be a unique probe for dark energy. 
On the other hand, such a NS-binary foreground should be subtracted to detect primordial GWs 
produced during inflation. Thus, the constraining power on dark energy and the detectability of 
the primordial gravitational waves strongly depend on the detector sensitivity and are in close re- 
lation with one another. In this paper, we investigate the constraints on the equation of state of 
dark energy with future space-based GW detectors with/ without identifying the redshifts of host 
galaxies. We also study the sensitivity to the primordial GWs, properly dealing with the residual 
of the NS-binary foreground. Based on the results, we discuss the detector sensitivity required to 
achieve the forementioned targeted study of cosmology. 

PACS numbers: 



I. INTRODUCTION 



Future space-based gravitational- wave (GW) detectors 
such as DECI-hertz Interferometer Gravitational-wave 
Observatory (DECIGO) d H and Big-Bang Observer 
(BBO) 3 (see also 3 for updated information) arc the 
most sensitive to GWs in 0.1 — 1 Hz band and will aim 
at detecting the primordial GW background, the merg- 
ers of intermediate- mass black holes (BH), and a large 
number of neutron-star (NS) binaries in an inspiraling 
phase. These GW sources enable us to probe inflation 
and density perturbation in the early universe [^-Q, to 
measure the cosmic expansion with unprecedented preci- 
sion ITo| , to investigate the population and formation 
history of compact binary objects and to test alter- 
native theories of gravity |12h15|. Therefore, DECIGO 
and BBO will open up the window of gravitational-wave 
cosmology. 

It is known that the continuous GW signal from a 
compact-binary object provides a unique way to mea- 
sure the luminosity distance to the source with high pre- 
cision. Such binary sources are often referred to as the 
standard siren (analogous to the electromagnetic stan- 
dard candle). With the redshift information determined 
by an electromagnetic follow-up observation, the stan- 
dard siren can be an accurate tracer of the cosmic ex- 
pansion [iHl- The potential power of this method as a 
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dark-energy probe has been investigated with ground- 
and space-based detector configurations for GWs [1, [TtI - 
[25} . Although these detectors will constrain the evolution 
of dark-energy, most of the preceding works assume that 
the redshifts of all GW sources are known by electro- 
magnetic observations. This assumption is rather strong 
and is too optimistic to be justified, because the spec- 
troscopic follow-up observation of galaxies is, in general, 
time-consuming, particularly at high redshifts 26[. In 
addition, not all galaxies can be observed due to intrin- 
sic faintncss, the absence of spectral features, limited sky 
coverage, and a limited redshift range (redshift desert). 
As a result, in a practical follow-up observation, the 
fraction of the binary sources whose redshifts are spec- 
troscopically obtained is significantly reduced. From a 
rough estimate based on the number density of galaxies 
potentially observable and the number of galaxies to be 
observed in the future galaxy redshift surveys, it turns 
out that the fraction of redshift identification could be 
~ 10^^ with uncertainty of about one order of magni- 
tude (See Sec. |VF]in detail). This means that the GW 
sources with redshift information would be rare unless 
we perform a large-scale follow-up campaign dedicated 
for GW events. Thus, this reduces the number of binary 
sources available as a standard siren. 

Another approach is to measure the cosmic- 
expansion rate from GW observations alone without 
electromagnetically-estimated redshifts. As suggested in 
the cosmic expansion affects not only the amplitude 
(luminosity distance) of GWs but also the phase, and it 
is possible to directly obtain information about the cos- 
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mic acceleration by accurately measuring the GW-phase 
shift of a binary source at a certain redshift. Although 
the redshifts of the binaries are assumed to be deter- 
mined by the electrornagnetic follow-up observations in 
the previous works [ll.l27l|. we show in this paper that we 
can measure the expansion history of the universe with- 
out any reference to the electromagnetic counterpart or 
host-galaxy identification by combining the luminosity 
distance and the cosmological phase shift, which are both 
independently determined from the amplitude and phase 
of GWs. This method based on a purely GW observation 
enables us to compare the observational data with those 
obtained in other electromagnetic observations. Since 
proposed space-based GW detectors such as DECIGO 
and BBO would detect ^ 10^ NS binaries, they provide 
a novel opportunity to measure the property of dark en- 
ergy without using the cosmic ladders [J| . 

Our primary interest here is in the potential of the 
standard siren without any help from redshift informa- 
tion of the binary sources. However, we also show how 
the sensitivity is improved by adding the redshift infor- 
mation of a fraction of NS binaries. As a result, we 
found that the contribution of binary sources whose red- 
shift is determined tightly constrain the dark-energy pa- 
rameters if the fraction of redshift identification is larger 
than 10-5 - lO"'^. In some cases that the fraction of 
redshift identification by the future galaxy redshift sur- 
vey is smaller than the one we expect, the contributions 
of binary sources with/ without redshifts arc compara- 
ble. Therefore, our study in this paper indicates that the 
GW standard siren without redshift information guaran- 
tees minimally achievable constraint on the dark energy 
and that the follow-up observation of the host galaxies 
targeted at the GW sources is crucial to improve the sen- 
sitivity. 

Another important scientific target of DECIGO and 
BBO is a primordial GW background. For its search, 
one needs to accurately identify the waveforms of NS 
binaries and subtract them up to the level sufficiently 
below the amplitude of the primordial GW background 
[28l [29j . Thus, an accurate determination of the wave- 
forms of NS binaries is essential both for the utilization 
of the NS binaries as the standard siren and the detec- 
tion of the primordial GW background. For this reason, 
in this paper, we comprehensively treat these issues and 
investigate the scientific outcome obtained by DECIGO 
and BBO, allowing their noise curves to be scaled ap- 
propriately and including confusion noises produced by 
astrophysical sources. The scaling of the noise curve is 
worth considering, because the requirement for instru- 
mental noise is fiexible at the preconceptual stage of the 
detector configuration and should be determined to op- 
timize the scientific results. 

This paper is organized as follows. In Sec. [Hi we 
briefly review the detector configurations of DECIGO 
and BBO and provide a model for the power spectra of 
detector noises and astrophysical confusion noises that 
we use throughout this paper. In Sec. IIIIl we present the 



model of a GW waveform from a compact binary and give 
an estimate of the signal-to-noise ratio obtained when a 
matched-filtering analysis is performed. In Sec. IIVI we 
explain the procedure of the NS-binary subtraction and 
give the residual contribution to the total noise after the 
subtraction. Based on these results, we calculate the de- 
tector sensitivities to dark energy in Sec. |V] and to a 
primordial GW background in Sec. IVII In Sec. I VIII we 
give discussions on the results and the feasibility of each 
method. Sec. I Villi is devoted to conclusions. Throughout 
the paper, we adopt units c = G = 1. 



II. DETECTOR NOISE AND ASTROPHYSICAL 
FOREGROUNDS 



In this section, we summarize the detector configura- 
tion of DECIGO and BBO and provide detector-noise 
curves that are used throughout the paper. 



A. DECIGO/BBO 

In the current preconceptual design of the detector 
lll, DECIGO and BBO orbit the Sun with a period of one 
sidereal year, and constitute four clusters, each of which 
consists of three spacecrafts exchanging laser beams with 
the others. Two of the four clusters are located at the 
same position to enhance the correlation and hence the 
sensitivity to a stochastic GW background, and the other 
two are widely separated from each other on the Earth 
orbit in order to enhance the angular resolution so that 
we can easily identify the host galaxy of each NS binary 
via the electromagnetic follow-up observations [1] . 

DECIGO and BBO arc the most sensitive to GWs in 
the 0.1 — 1 Hz band, which is determined so as to avoid the 
astrophysical GW foreground produced by white-dwarf 
(WD) binaries. DECIGO instrumental noise is domi- 
nated by radiation pressure noise below ^ 0.2 Hz and 
by laser shot noise above that frequency. As for BBO, 
the sensitivity at high frequencies is limited by beam- 
pointing-jitter and stray light noises rather than laser 
shot noise. The noise curve of DECIGO is calculated 
adopting currently-proposed design parameters 0, [30| . 
We obtained fitting formula for the sky-averaged noise 
curve of DECIGO single interferometer: 



^rD(/) = 3.30x10- 



-1-3.09 X 10" 



-50 



(iHz) 



Hz-i , (1) 



where /c = 7.69 Hz. On the other hand, BBO sky- 
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averaged noise curve is given in Q as 

SrB(/)= 6.15 X 10-51 (^^^ +1.95x10- 

+ 1.20 X 10"''*^ ( Hz-i . 
V IHz/ 



(2) 



DECIGO and BBO have different noise shapes, result- 
ing from the difference in the interferometer type and op- 
tical parameters. The default noise curve of DECIGO is 
designed to avoid WD confusion noise at low frequencies. 
On the other hand, BBO has slightly better sensitivity 
than DECIGO in the frequency range below several Hz, 
but the WD confusion noise dominates the noise curve 
below ~ 0.1 Hz. As we will quantitatively discuss later 
including astrophysical confusion noise, the BBO sensi- 
tivity in amplitude to a NS binary is nearly three times 
better than that of DECIGO. 



B. Astrophysical foregrounds 

DECIGO and BBO are sensitive to a large number of 
astrophysical GW sources, in particular NS, BH, and WD 
binaries. It is expected that the low-frequency side of the 
noise curve would be dominated by the GWs from cosmo- 
logical (extra-galactic) population of WD binaries, which 
still remains after subtracting individually identified sig- 
nals in the process of data analysis and are stochastic in 
nature. According to the estimation by |3l| . the fitting 
formula for the power spectrum of WD confusion noise 
[3^ (residual contribution after the subtraction) is given 

by 



5WD(/)^4.2x 10- 



X exp 



IHz 
/ 



-7/3 



5 X 10-2 Hz 



Hz- 



There is also contribution from galactic WD binaries. 
However, as pointed out in |33j . the galactic sources are 
sufficiently sparse in frequency space above 3 x lO^"^ Hz 
so that one is able to fit them out of the data. Actu- 
ally, using the expression for the residual foreground after 
the subtraction given in [3^ . we verified that the galac- 
tic sources contribute less than 10 % to the total-noise 
curve above the frequency 5 x 10~^Hz. Thus, galactic 
contribution can be ignored for our present purpose. 

Another astrophysical source that we have to take into 
account is the NS-binary foreground. According to jssj . 
the energy density of GWs from NS binaries per loga- 
rithmic frequency bin normalized by the critical energy 
density of the universe at present is written as 



9H^ 



no = 



n{z) 



(1 + z)4/3ii-(z 



■dz 



(3) 



Hq is the Hubble constant and Mc is the chirp mass de- 
fined as Mc = if/^Mt, together with the total mass Mt = 
mi 4- m2 and the symmetric mass ratio ?/ ~ mim2/Mf . 
'h(z) is the NS merger rate per unit comoving volume per 
unit proper time at a redshift z. We adopt the following 
fitting form of the NS-NS merger rate given in [2^ : 



n{z) = ho s{z) ; s{z) = { |(5-z) (1< z < 5) 




(4) 

where the function s{z) is estimated based on the star 
formation history inferred from the UV luminosity [s^ ]. 
The quantity no represents the merger rate at present. 
We assume the flat ACDM universe with il^ = 0.3 = 
0.7) and the Hubble parameter H{z) is given by 



H{z) = Ho [n^{i + zf + (1 - n^)] 



1/2 



(5) 



Writing Ho = h-j2 x 72 km s ^ Mpc ^ and substituting 
Eqs. dH) and jS]) for Eq. dS]), we obtain 



^l^^^U) = 3.74 X 10 ,.72 



Mr 



1.22M, 



5/3 



IHz 



2/3 



no 



10-^Mpc-''yr 



which is translated into the NS confusion-noise power 
spectrum 



sru) = 1.55 X 10-4^2' 



Mr \ V / ^ -^'^ 



,1.221/0/ VlHz 
using the relation 



no 



10-6Mpc-Vi-i 



3Hl 
47r2/ 



3 s 



(6) 



BH-BH and NS-BH binaries are also the sources of the 
foreground. However, the populations are smaller, com- 
pared with NS-NS binaries. Therefore we neglect them 
for simplicity. 



C. Total noise curve 

DECIGO and BBO have almost the same sensitivities, 
but BBO is approximately three-times more sensitive to a 
GW from a NS binary than DECIGO in terms of signal- 
to-noise ratio (SNR) (e.g. for a NS binary at the red- 
shift z = 5, sky-position- and inclination-averaged SNR 
is 8.5 and 26.3 for a single interferometer of DECIGO 
and BBO, respectively.). In other words, the DECIGO 
noise curve improved by a factor 3 in amplitude is com- 
parable to that of BBO. Our purpose in this paper is to 
assess cosmology achieved by space-based detectors such 
as DECIGO and BBO, and to make the requirement for 
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FIG. 1: Noise curves of default DECIGO (solid blue curve) 
and those scaled by r„ = 1/2 (green), 1/3 (orange), 1/5 (red), 
which include the confusion noise from a number of WD bi- 
naries. The WD-binary foreground \/S^^ is shown with 
solid, black curve. T he th ree diagonal lines represent the NS- 
binary foreground \/S^ before subtraction (or 7?.ns ~ 1) 
with no = 10"^ 10"", 10"'^ Mpc"^ yr"^ from the top to the 
bottom. 



the experimental design clear. Thus, we consider only 
DECIGO, but allowing the noise curve in Eq. ^ to vary 
by an overall factor r„ = 1, 1/2, 1/3, 1/5 in amplitude. 
Namely, the power spectrum scales as 

Note again that BBO noise power spectrum in Eq. ([2]) 
approximately corresponds to that of DECIGO with r„ = 
1/3. 

Including astrophysical contributions, a total-noise 
power spectrum is given by 



NS 



(7) 



Here the factor T^ns denotes a suppression factor due 
to the subtraction of individually identified NS binaries, 
which will be estimated in Sec. IIVI The contribution of 
each term in Eq. ([7]) is shown in Fig. [TJ At high frequen- 
cies, roughly above 0.2 Hz, the noise curves scale straight- 
forward with r„ if we neglect the contribution from NS 
binaries. On the other hand, at low frequencies, the WD 
confusion noise prevents the total noise from improving 
as r„ decreases. This means that, from the experimental 
point of view, improvement required for the instrument is 
involved with not radiation-pressure noise but laser shot 
noise or beam-pointing-jitter and stray light noises. 



and primordial GW backgrounds. Firstly, in Sec. IIVI we 
compute SNR for each NS binary with the noise power 
spectrum, S^^^^{f) + S^^{f), assuming that large sig- 
nals well above the noise curve are easily subtracted. 
Then, we determine how much the residual amplitude 
of the NS-binary foreground remains after the subtrac- 
tion process, which depends on the SNR threshold value. 
In the next step, the noise spectrum, S'il'''\f ) + S^^{f) + 
S'^^(/)72,NSi is used to calculate the Fisher matrices 
for parameter estimation of NS binaries and cosmolog- 
ical models. The contribution of the confusion noise, 
S^^if) + S^^{f)'JlNs, is also taken into account in the 
correlation analysis for a stochastic GW background. 



III. GW WAVEFORM OF A BINARY AND SNR 

For a single binary system, the Fourier transform of the 
GW waveform is expressed as a function of frequency / 



Hf) 



A 



dUz) 



(8) 



where is the luminosity distance, and the quantity 
AIz = (1 -I- z)AIc is the redshifted chirp mass. Ale is 
the proper chirp mass defined in the source rest frame. 
The constant A is given hy A ~ (a/6 tt^/'^)"^, which in- 
cludes the factor •\/4/5 for a geometrical average over 
the inclination angle of a binary [s^]. The function 4'(/) 
represents the frequency-dependent phase arising from 
the orbital evolution, and at the order of the restricted 
1.5 post-Newtonian (PN) approximation, it is given by 

mm 



25 



(9) 



where rj = mim2/ (rni + m-i)^ is the symmetric mass 
ratio, and ic and are the time and phase at coales- 
cence, respectively. The first term in the square brackets 
in Eq. ^ corresponds to the Newtonian-order dynam- 
ics and the second and third terms represent the post 
Newtonian-order corrections in powers of t; '--^ (ttM^/)^/'^. 
The last term in the bracket is a phase correction due to 
cosmic expansion [l|,|23, where X{z) is defined as 



D. Outline of analysis and noise model 



X(z).iU-^ 
^ ' 2 I " 1 + z 



In this paper, we deal with the instrumental noise and 
the astrophysical confusion noise in a self-consistent man- 
ner when we estimate the sensitivities to dark energy 



or equivalently expressed as [a(0) — d(z)]/2. Thus, 
X{z) > {X{z) < 0) corresponds to the accelerating 
(decelerating) universe. Note that this correction is "—4 
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PN" order and becomes more important at lower fre- 
quencies. This is because a binary longer stays at low 
frequencies (when the binary separation is larger) and 
has the larger cycle number. 

The squared SNR of a binary GW signal at a redshift 
z is defined by 



,=1 



(10) 



Here 5'h(/) and /i(/) are given by Eqs. ([7]) and dH), re- 
spectively. The summation in Eq. pOp is taken with 
respect to the number of independent detectors of DE- 
CIGO. (DECIGO has four clusters, each of which has 
two independent interferometers). Since it is assumed 
that all interferometers have identical noise-levels, the 
summation can be replaced just with an overall factor of 
8. The lower cutoff frequency /min is given by the func- 
tion of observation time Tobs as well as the redshift and 



/n 



0.233 



"m7 



5/8 



To 



bs 



3/8 



Hz 



(11) 



The upper cutoff of the frequency /max naturally arises 
from the noise curve, since the coalescence frequency of 
the binary with the mass Mq is typically ~ kHz. For 
the computational purpose, we set /max = 100 Hz. 



IV. NEUTRON-STAR BINARY SUBTRACTION 

NS-binary foreground is not stochastic in the DECIGO 
frequency band [28| and in principle can be subtracted 
from observation data. How well one can subtract them 
and how much the residuals are left depend on the esti- 
mation accuracy of binary parameters, or the SNR of the 
signal from each NS binary. This subtraction process is 
mandatory to use NS binaries as a probe for dark energy 
and to search for a primordial GW background. 

In a recent paper [l^, Yagi and Seto have performed 
Monte-Carlo simulations generating synthetic GW data 
for binaries at random sky positions and with orienta- 
tions of angular-momentum axes over redshifts up to 
z = h. They estimated how much residuals are left after 
the subtraction process. The criterion for the successful 
subtraction is that SNR pi for each binary exceeds the 
SNR threshold pth- The fraction of the residuals in the 
energy density of NS binaries after the subtraction can 
be expressed as 



F[z)^ 
together with 



l^dzF{z)mll{!^z) 
l^dzl^^l^{],z) 

E»p|(^)Qbth-p»(g)] 



5^,,/-5/3 .2/3 

ml ^ ^ 



h{z) 




0.6 



0.8 



1.0 1.2 1.4 

Ah/ A (5) 



1.6 



{i + zy/^H{z) 



FIG. 2: Fraction of NS residual 7?,ns after the subtraction as 
a function of pth/pe(5). The fitting formula is presented in 
Appendix 1X1 



where 0[- • •] is the step function. 

In Ref . [231 , the fractional residual 7?,ns is provided as 
a function of pth/Pc(5), where Po(5) is the SNR for an 
edge-on NS binary at the redshift z ~ 5. The 7?.ns is 
shown in Fig. ^ The figure indicates that T^ns rapidly 
decreases at the point where pth/Pe(5) ~ 0.79 or cquiva- 
lently Po(5) is nearly 30% better than pth- In Appendix 
El based on this result, we give an explicit fitting form 
of the result for the later calculation of 7?.ns • 

The pe(5) is related with the inclination-angle averaged 
SNR p(5) defined in Eq. HU]) by po(5) = %/5/4 x p(5). 
This is because the angular dependences of GW wave- 
forms are /i+ cx (1 -|- cos^6')/2 and oc cos 9 for each 
polarization mode, where is the angle between the line 
of sight and the axis of the orbital motion of the bi- 
nary. In order to calculate p(5), we use the noise power 
spectrum in Eq. ([7]) without the NS-binary contribution. 
Strictly speaking, the NS-binary foreground should be 
initially included and each binary signal has to be iter- 
atively subtracted one by one for the consistency of the 
calculation. However, as the authors of [2^ have per- 
formed, we can simplify the procedure by neglecting the 
contribution from the NS-binary foreground in the calcu- 
lation of SNR. This simplification would be valid because 
large GW signals can be easily subtracted from the data 
to the noise level that is limited by the instrumental noise 
and the WD foreground. We choose the SNR threshold 
as pth = 20. 

Once p(5) is obtained, we can calculate pth/Po(5) 
and translate it into 72.ns by using the fitting for- 
mula in Appendix [X] The values of 72,ns for observa- 
tion time, Tobs = 3, 5, 10 yr, and noise-curve scaling, 
r„ = 1, 1/2, 1/3, 1/5, are listed in Table H] If r„ is less 
than a half, all NS binaries are fitted out. Even in the 
case of r„ = 1, t he am plitude of the residual of the NS 
foreground is ^ V2OO ~ 14 times smaller than the initial 
foreground levels and hardly affect the shape of the noise 
curves presented in Fig. [T] Note, however, that the sup- 
pression of the NS foreground below the noise curve is not 
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always sufficient for tlic detection of a GW background, 
because the residual may contaminate a correlation sig- 
nal. This issue is discussed in Sec. IVll 



V. STANDARD SIREN AS A PROBE FOR 
DARK ENERGY 

A. Standard siren 

The continuous GW signal from a compact-binary ob- 
ject provides a unique way to measure the luminosity dis- 
tance to the source with high precision [l6j . Such binary 
sources arc often referred to as standard sirens. There 
are two different aspects in using the standard siren. 

One is to measure the amplitude of GWs, or the lu- 
minosity distance as a function of redshifts (see Eq. ([8|)). 
However, GW information alone cannot provide a source 
redshift, since the mass parameter we can determine is 
Mz = {l-\-z)Mc and the source redshift degenerates with 
the proper source mass Mc [1^. So the redshift has to 
be determined from electromagnetic observation of the 
host galaxy. This method as a tool to measure the cos- 
mic expansion has been investigated by many authors 
0j [l3423l J assuming that source redshifts are known by 
spectroscopic follow-up observations. Thus, the feasibil- 
ity of using NS binaries as the standard siren relies on 
the determination of the redshift of each binary, which 
requires high-angular resolution for GW detectors to se- 
lect a true host galaxy out of several candidates. A lack 
of angular resolution of ground-based detectors J3£ 
or Laser Interferometer Space Antenna (LISA) [li 
limits the availability of the standard sirens as a cosmo- 
logical probe, though the identification of the host galaxy 
is possible for the binary event accompanying an elec- 
tromagnetic counterpart. On the other hand, accord- 
ing to Cutler and Holz BBO has the angular reso- 
lution ^ 1 — lOOarcsec^ and can uniquely identify the 
host galaxy of the binary. The authors also have shown 
that cosmological parameters can be measured with high 
precision, using a large number of NS binaries, ^ 10^. 

Another method utilizes the phase modulation of GWs 
due to the cosmic expansion, which is the last term in- 
cluding X{z) in Eq. ® P,[23. Since the amplitude and 
phase bring us independent information, this method is 
complementary to the method described above. While 
the redshifts of the binaries are assumed to be determined 
by the electromagnetic follow-ups in the preceding works 
Ul, : we point out in this paper for the first time that 
we can measure a cosmic-expansion rate only with GW 
observations without electromagnetically-estimated red- 
shifts by combining the luminosity distance and the phase 
shift. This is a great advantage of our method, because a 
pure GW observation enables us to compare its observa- 
tional data with those obtained in other electromagnetic 
observations. However, unfortunately, this method is not 
applicable to advanced ground-based detectors [42| - |45j 
and LISA as a sensitive cosmological probe, because of 




the weakness of the signal coming from the cosmic expan- 
sion. To detect the signal, a large number of GW sources 
is necessary. DECIGO and BBO are expected to detect a 
lot of NS binaries, which enable them to detect the weak 
signal. In the subsequent sections, we first investigate the 
power of utilizing the phase correction, without assuming 
any electromagnetically-estimated redshift, and give the 
figure of merits. Then we will compare this method with 
the previous one, in which source redshifts are assumed 
to be given, and combine both methods to maximize the 
accuracy of the estimation. 



B. Constraints on cosmological parameters 

For simplicity, we consider a spatially flat universe with 
dark energy whose equation of state is parametrized as 
w{z) ~ Wf) + Waz/{1 + z) [4^. Then the luminosity dis- 
tance and the cosmological phase correction in Eqs. 
and (13 are written as 



dL{z) = (1-1-2:) 



dz' 
W)' 

H{z)^Ho{nrail + zf 

+ {1 - n„,)(l + z)3(i+"'o+'^") exp 



SWnZ 



1 + z 



(12) 
(13) 

1/2 



In the above, wq, Wa, i^m, and Hg are the parameters 
to be determined from GW observations. For the anal- 
ysis of error estimation, we adopt a fiducial set of cos- 
mological parameters: wq = —1, Wa = 0, flm = 0.3, 
Hq = 72kms^^ Mpc~^. For illustrative purpose, we give 
some dL - X plots for certain parameter sets in Fig. |31 
The redshift dependence of the curves are different de- 
pending the values of wq and Wa so that in principle one 
can independently determine these parameters. 

In what follows, we calculate the estimation errors of 
binary parameters 9a'- Mz, t], tc, ipc dL, and X in the 
waveform of Eq. ([S]). Then we derive the estimation er- 
rors of the cosmological parameters from them. The fun- 
damental tool wc use to estimate the errors for a single 
binary is the Fisher matrix formalism. The Fisher matrix 
for a single binary is given by [s^ \^ 



/n 



dah*if)dbhif) 



df 



(14) 



where da denotes a derivative with respect to a parame- 
ter 9 a- Since the eight interferometers of DECIGO and 
BBO are identical, the summation is reduced to just mul- 
tiplying by a factor of 8. The noise power spectrum and 
GW signal are given by Eqs. ([7]) and ([U, respectively. 
The frequency cutoffs are also given by Eq. (|lip and 
/max — 100 Hz. Given the numerically evaluated Fisher 
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0.0284 Hz) 




P(5) 


Pth/Pe(5) 




P(5) 


Pth/Pe(5) 


7?.NS 


P(5) 


Pth/Pe{5) 


TZns 


1 


24.06 


1.487 


4.66 X 10"^ 


24.07 


1.486 


4.64 X 10"^ 


24.09 


1.485 


4.62 X 10"^ 


1/2 


47.75 


0.749 





47.77 


0.749 





47.78 


0.749 





1/3 


71.25 


0.502 





71.25 


0.502 





71.26 


0.502 





1/5 


118.9 


0.304 





117.9 


0.303 





117.9 


0.303 






0.10 



0.05 



0.00 



-0.05 



-0.10 



-0.15 



-0.20 



TABLE I: Fraction of NS residual 7?.ns for each parameter set of Tobs and r„. 
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-1 


Wo = 


-0.8 
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10 



-0.10 



-0.15 



-0.20 




Wa 


= 


Wa 


= 0.5 


Wa 


= 1.0 



10 



dL{z)IH^ 



-1 



FIG. 3: dL - X plot as a function of redshift. Left: wa is varied, while Wa and Sim are fixed. Each curve corresponds to 
wo — —1.2 (red, dotted curve), wo — —1 (green, solid curve), wo ~ —0.8 (blue, dashed curve). Right: Wa is varied, while wo 
and Qn-i are fixed. Each curve corresponds to Wa = (green, solid curve), Wa ~ 0.5 (red, dotted curve), Wa = 1.0 (blue, dashed 
curve). The green solid curves in both figures are for ACDM universe. The points corresponding to z — 1,2,3,4,5 are shown 
on the curve. 



matrix, the marginalized 1-a error of a parameter, A9a, 
is estimated from the inverse Fisher matrix 



by |2g 



A9a = v/{r-l}aa. 

Once the estimation error of X is obtained, it is 
straightforward to calculate the estimation errors of the 
cosmological parameters. As described in Appendix |b1 
the error of is much smaller than that of X. This 
fact means that we can replace an observed cIl with the 
corresponding redshift in the fiducial cosmological model 
when we derive the measurement accuracies of cosmolog- 
ical parameters from X. Thus, for the simplicity of the 
analysis, we use X{z) instead of Ar(di). Futhermore, we 
assume that the Hubble constant Hq is known a priori, 
and fix it to Hq — 72 kms~^ Mpc~^, because the Hubble 
constant has been determined at a-few-percent level from 
the observation of nearby Ccphcids and supernovae j48| . 
Thus, the free parameters of the Fisher matrix are wo, 
Wa, and fljri- The Fisher matrix is given by 



Fab 



daX{z)dbX{z) dN{z) 



dz 



dz , 



(15) 



where dN{z)/dz is the number of NS binaries in the red- 
shift interval [z, z + dz] observed during Tobs and is given 



dN{z) 
dz 



47rr^(z) h{z) 
H{z) 1 + 2 



(16) 



where r[z) is the comoving radial distance defined as 
r{z) = dL{z)/(\ + z) and fi is given in Eq. Since 
the normalization of h is still uncertain, we adopt the 
most recent estimate, no = 10~^ Mpc~'^ yr~^, as a reli- 
able estimate based on extrapolations from the observed 
binary pulsars in our Galaxy [49j . and also consider op- 
timistic and pessimistic values, 10^^, lO^*" Mpc""^ yr^^. 
Since the observation time Tobs is also a crucial parame- 
ter, we suppose 3-yr, 5-yr, and 10-yr observation for the 
Fisher matrix. In Eqs. ((TS|) and ((T5)) . the integration is 
performed explicitly with respect to a redshift. Note, 
however, that the integrand can be regarded as a func- 
tion of the luminosity distances, instead of the rcdshifts, 
as discussed in Appendix |B] 



C. Results 

Since we are interested in parameters involved with 
dark energy, the measurement accuracies of wq and Wa 
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Tohs = 3 yr 




no = 10"' 


no = 10"" 


no = 10"" 


rn 


Awo 


AWa 


FoM 




AWa 


FoM 


Atoo 


AWa 


FoM 


1 

1/2 
1/3 
1/5 


3.05 X 10" 
1.53 X 10° 
1.03 X 10° 
6.26 X 10-1 


2.63 X 10' 
1.37 X 10^ 
9.55 X 10" 
6.15 X 10" 


1.71 X 10"^ 
6.32 X 10"^ 
1.31 X 10"' 
3.20 X 10"' 


9.81 X 10"' 
5.21 X 10"^ 
3.25 X 10"' 
1.98 X 10"' 


8.44 X 10" 
4.59 X 10° 
3.02 X 10" 
1.95 X 10" 


1.66 X 10"' 
5.65 X 10-1 
1.31 X 10" 
3.20 X 10" 


3.59 X 10"' 
1.53 X 10"' 
1.03 X 10"' 
6.26 X 10-^ 


2.99 X 10° 
1.37 X 10° 
9.56 X 10"' 
6.15 X 10-' 


1.33 X 10" 
6.32 X 10" 
1.31 X 10' 
3.20 X 10' 



Tohs = 5 yr 




no = 10"' 


no = 10"" 


no = 10"" 


r„ 


Awo 


AWa 


FoM 


AuiQ 


AWa 


FoM 


Awo 


AWa 


FoM 


1 

1/2 
1/3 
1/5 


1.03 X 10° 
5.24 X 10"' 
3.59 X 10-' 
2.28 X 10"' 


1.02 X 10' 
5.46 X 10° 
3.89 X 10° 
2.61 X 10° 


1.18 X 10-' 
4.10 X 10"' 
8.12 X 10-1 
1.82 X 10° 


3.30 X 10-' 
1.66 X 10"' 
1.13 X 10-' 
7.20 X 10-^ 


3.25 X 10" 
1.73 X 10" 
1.23 X 10° 
8.26 X 10-' 


1.14 X 10" 
4.10 X 10" 
8.13 X 10° 
1.82 X 10' 


1.20 X 10"' 
5.24 X 10"^ 
3.59 X 10"^ 
2.28 X 10"^ 


1.14 X 10" 
5.46 X 10-' 
3.89 X 10-1 
2.61 X 10"' 


9.21 X 10" 
4.10 X 10' 
8.12 X 10' 
1.82 X 10^ 




Tobs = 10 yr 




no = 10"' 


no = 10"" 


no = 10"" 




Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


1 

1/2 
1/3 
1/5 


2.57 X 10"' 
1.40 X 10-' 
1.03 X 10"' 
7.28 X 10-^ 


2.84 X 10° 
1.64 X 10° 
1.23 X 10° 
8.99 X 10-' 


1.56 X 10" 
4.73 X 10" 
8.33 X 10" 

1.57 X 10' 


8.25 X 10"^ 
4.43 X 10-^ 
3.24 X 10"^ 
2.30 X 10"^ 


9.06 X 10"' 
5.17 X 10-' 
3.90 X 10-' 
2.84 X 10-' 


1.53 X 10' 
4.72 X 10' 
8.33 X 10' 
1.57 X 10^ 


2.94 X 10"^ 
1.40 X 10-^ 
1.03 X 10"^ 
7.28 X 10-^ 


3.15 X 10"' 
1.64 X 10-' 
1.23 X 10-' 
8.99 X 10-^ 


1.26 X 10^ 
4.73 X 10^ 
8.33 X 10^ 
1.57 X 10^ 



TABLE II: Measurement accuracy of wo and Wa, and FoM defined in Eq. (|18fl for each parameter set of Tobs, n-o, and r„. In 
the table, no is in the unit of Mpc"'' yr"'. 

in the phase-correction term. So, from Eq. (jl4p . 
a"^ - r(--'^) oc / iLUlL^^^Hdf 

.-23/3 rT,23/8 
^ Jmin ^ ^obs ' 

In the second Une of the above equation, we used 
/mill oc T^T^J^. In total. Tab oc xfjj/*, which leads to 

Awo, Aioq cx Tobs^^^^- l^otc that the r„ dependence in 
Eq. (jl7p is not a good approximation for large Tobs be- 
cause of the presence of the WD foreground. For the 
same reason, the deviation of the above scaling for Tobs 
is also large for small r„. Thus, the scaling relations in 
Aiuo and Awa should be considered as qualitative one. 

To show the results more clearly, let us define the figure 
of merit (FoM) for dark energy. 



FoM = Vdet -iab , (18) 

where ^ab is the inverse matrix of (r-^)^;, with a,b ~ 
wq, Wa- As is obvious from the definition, FoM is in- 
versely proportional to the area of an error ellipse on the 
Wo - Wa plane. The FoM is listed in Table HIl and is plot- 
ted in Fig. [5] as a function of r„ for various values of ho 
and Tobs- From the scaling of Awq and Awa, FoM is 
expected to scale as FoM cx no Tf^J/^r"^. 



marginalized over the other remaining parameter are tab- 
ulated in Table HIl They change significantly depending 
on the parameters, no, Tobs, and r„. For example, if we 
choose the intermediate values, no = 10"® Mpc"'^ yr"i 
and Tobs = 5 yr, the errors Awo and Awa are 0.330 and 
3.254 for r„ = 1 and 0.113 and 1.231 for r„ = 1/3. For 
10-yr observation, Auio and Awa reach at the levels of 
several xlO-^ and a few xlO-^, respectively. Error el- 
lipses on the Wo - Wa plane are shown in Fig. |4] for the 
fixed no varying the values of r„ and Tobs- 

Scaling relations of the measurement accuracy with no, 
Tobs, and r„ help us understand the results. Awq and 
Awa are expected to scale like 

Awq, Awa oc Vt^ (X Tig ^^%bf^^^r„ . (17) 

This relation can be understood as follows. The no de- 
pendence purely comes from the number of NS binaries 
through h{z) in Eq. (|16p . The r„ dependence is in- 
ferred from the noise spectrum in Eq. (|14p . which is 
proportional to r^. The observational-time dependence 
is rather complicated and comes from two factors: the 
number of NS binaries and the GW phase correction in- 
volved with X. From Eq. Tab oc dN/dz oc Tobs- 
On the other hand, from Eq. the cosmological- 
phase-correction term is larger at low frequencies and 
is roughly proportional to f^l^^^ ■ Around the frequency 
/mill, noise is dominated by the WD foreground, which is 
Sj^^ oc /-^/'^ and effectively reduces the power of /min 
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FIG. 4: Wo - Wa error ellipses marginalized over Q,m- The merger rate is fixed to ho = 10"*^ Mpc""^ yr~^. From the left to the 
right panels, the observation time is Syr, 5yr, and 10 yr. In each panel, the larger to smaller ellipses denote those with r„ — 1, 
1/3, and 1/5, respectively. The dot at the center is our fiducial value: wo = —1 and Wa = 0. 
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FIG. 5: FoM achievable for each parameter set of r„, Tobs, and no. From the left to the right panels, the observation time 
is Syr, Syr, and 10 yr. In each panel, the merger rate is represented by curve types: ho = 10^^ (dashed), 10~^ (solid), 10~^ 
(dotted), in the unit of Mpc~'^ yr"^- 



D. Redshift dependence of measurement accuracy 

The measurement accuracies in the previous subsec- 
tion are calculated with NS binaries detected up to z = 5. 
From an observational point of view, it is interesting to 
see which redshift sources largely contribute to the mea- 
surement accuracies. To do this, we fix the parameters 
to Tn = 1/3, Tobs = 5yr, and no = IQ-^ Mpc'^ yr"! 
and compute the measurement errors and FoMs, restrict- 
ing the redshift range of the binary distribution below 
^max = 0.5 X i (i = 1, 2, • • • , 10), i.e. upper cutoff of the 
integral in Eq. (fTS)) . The results are shown in Fig. |6]in 
terms of the sensitivity fractions: FoM/FoM(Zfnax = 5), 
Awo(zmax = 5)/Awo, and Awa(zmax = 5)/Awa as a 
function of Zmax, where say, Au'o(2;inax = 5) represents 
the wq error with full redshift range. 

The figure shows that the sensitivity fractions sharply 
increase up to z = 1.5 — 2.0 and slowly approach the max- 
imum value at z = 5. This reflects merely the redshift 
distribution of NS binaries in Eq. In other words, the 



brighter source at low z gives large SNR but the number 
of the source is small. A large number of sources at mod- 
est z largely contributes the measurement accuracies of 
cosmological parameters. The interesting feature is that 
the slope of Aw;Q(zniax = 5)/Awa is gradual, compared 
with Ait;o(zmax = 5)/Ai(Jo, because the samples at high 
z are required to constrain the time variation of dark en- 
ergy, Wa (see Fig. El). The fraction FoM/FoM(2:max = 5) 
traces the redshift dependence of the Wa error and shows 
gradual increase. 

E. Adding redshift information of binaries 

In the previous subsections, we used the combination 
of the luminosity distance and the cosmological phase 
correction X to measure the dark-energy parameters. 
With this method, the cosmological expansion can be 
probed by GW observation alone without the identifi- 
cation of binary redshifts. Although it would be too 
ideal to assume that all binary redshifts are well deter- 
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FIG. 6: Redshift dependence of the measurement accuracies 
in d_L - X measurement: FoM/FoM(2;max = 5) (blue, solid), 

Awo(Zmax = 5)/Awo (red, dotted), A«;a(2max = 5)/AWa 

(green, dashed). Parameters are fixed to r„ = 1/3, Toba ~ 
Syr, and no = 10"** Mpc'^ yr~\ 



mined by electromagnetic follow-up observations, some 
fraction of them will be determined by the follow-up ob- 
servations. Then, we can optimize the constraints on the 
dark-energy parameters by using thus-determined red- 
shifts. In this subsection, we combine z - cIl and - 
X information and evaluate what fraction is needed to 
achieve FoM =100, which is a typical value expected in 
future dark energy surveys [sot . 

The relevant Fisher matrix is defined with the redshift- 
idcntification fraction a by 



dadL{z)dbdL{z) dN{z) 

daX{z)dbX{z) dN{z)^ 



dz 



'X 



(^) 



dz 



(19) 



where is an error in the luminosity distance, given 
in Eq. (jBl[) in Appendix |B] The free parameters of the 
Fisher marix are again wq, Wa, and fini, while the Hub- 
ble constant is fixed. Strictly speaking, the variables 
and X cannot be treated independently as in Eq. ([T^ . 
However, the correlation coefficient between and X 
is of the order of 10~^. This reflects the fact that am- 
plitude and phase are independent quantities, hence we 
can safely use the expression in Eq. ([T9|) . The fraction 
factor a would be a function of a redshift in a real galaxy 
survey, but we take it as a constant for simplicity of the 
analysis. 

In Appendix |B1 we estimate the error size of the lumi- 
nosity distance via the same procedure as in Sec. lVBl but 
including possible systematic errors: weak-lensing magni- 
fication due to the matter inhomogeneities along the line 
of sight and the peculiar velocity of each binary source. 
Then substituting the estimated errors for Eq. ([T9)) . we 
calculate a needed to achieve FoM =100, which we de- 
note by aioo, and the corresponding number of NS bina- 
ries Nz whose redshifts arc to be identified. 



10° 

o 

10^ 



^ Tohs = 5 yr 




no = iU 


















FoM=100 







FIG. 7: FoM as a function of a defined in Eq. or N^. 
The observation time and the merger rate are fixed to Tobs = 
Syr and no = 10~® Mpc"'' yr~^. The horizontal dashed line 
represents FoM = 100. The solid and dotted curves are the 
case with r„ = 1 and 1/3, respectively. 



Before showing the results, it would be helpful to see 
how FoM changes as a increases. In Fig. [3 FoM as a 
function of a is presented for fixed Tobs and ng. At low a, 
of course, the z - dL information plays no significant role 
in improving the FoM. But at a sa a few xlO~^, z - dL 
information begins to contribute to FoM, and enables the 
FoM to reach 100 at around a « a few xlO^^. The frac- 
tion aioo required to achieve FoM = 100 is listed in Ta- 
ble |TTT1 As anticipated, z - d^ measurement is much more 
powerful than the dL - X measurement because FoM|q=o 
is much smaller than FoM|c=i. However, how much frac- 
tion a is feasible in the future galaxy redshift surveys? In 
the next subsection, we estimate a from rough consider- 
ation based on the number density of galaxies potentially 
observable and the number of galaxies to be observed in 
the future galaxy redshift surveys. 



F. On the redshift determination from optical 
follow-up observations 

To obtain the redshift of a binary source, we need to 
identify its host galaxy and determine the redshift by 
an electromagnetic follow-up observation. As shown by 
Cutler and Holz [3] from the simple consideration based 
on the average number density of galaxies, DECIGO and 
BBO has the angular resolution ~ 1 — lOOarcsec^ and 
can uniquely identify the host galaxy of the binary. Al- 
though more complication such as galaxy clustering and 
so on should be taken into account, in what follows, we 
assume that a single host galaxy is found in a detector 
error cube. So our main concern here is what fraction of 
GW sources can we determine the redshifts by electro- 
magnetic follow-up observations in a realistic situation 
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70 
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3.0 X 10"^ 
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1.8 X 10"" 
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1.4 X 10"^ 
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Tobs = 10 yr 
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Qioo 




Qioo 




Qioo 




1 


2.6 X 10"" 


260 


2.2 X 10"" 


220 








1/2 


1.2 X lO""* 


120 


6.4 X 10"® 


60 








1/3 


9.4 X lO"'^ 


90 


1.5 X 10"® 


20 








1/5 


7.4 X lO"'^ 
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TABLE III: The fraction a required to reach FoM = 100 and the corresponding number of a NS binary that needs redshift 
identification. Tlie binary number is rounded off to tens, aioo = or A'^z = in the table means that FoM = 100 is already 
achieved by X measurement without any redshift information, no is in the unit of Mpc"^ yr"^. 



with finite observation time and limited survey magni- 
tude. 

The number density of galaxies potentially observable 
(brighter than a limiting apparent magnitude) can be es- 
timated from the galaxy luminosity function. According 
to the paper of the Hubble Ultra Deep Field [5l|, the 
total number of galaxies potentially observable per unit 
square angle with the limiting magnitude mum ~ 29 in 
an optical band is ~ 5 x 10^ /arcmin^. The total number 
over the celestial sphere is ~ 7x 10^^. However, in a real 
observation, various factors such as assigned observation 
time and accessible observation bands restrict the num- 
ber of galaxies to be observed in the survey. The future 
galaxy spectroscopic surveys, JDEM/WFIRST [s^ and 
Euclid [53, plan to observe ^ 10^ galaxies in the redshift 
range 0.5 < z < 2 with the redshift precision better than 
0.1% d^. Hence, the fraction of the galaxies whose red- 
shifts are listed in a galaxy catalog is /catalog 1 x 10"''. 
For simplicity, we assume that GW events randomly oc- 
cur in one of the galaxies. Unless we arrange follow-up 
observations especially dedicated for the GW events, the 
probability that the host galaxy of a GW event is listed 
in the galaxy catalog is also /catalog, which gives the frac- 
tion of redshift determination of binary sources. 

The estimation based on the galaxy catalog is sim- 
ple but seems to be reasonable, because galaxy redshift 
surveys are efficient in that they first select the candi- 
date galaxies for the spectroscopic observation by lumi- 
nosity and emission lines. This means that other galax- 
ies not listed in the catalog is more difficult to obtain 



the redshift. So we expect that the number of redshift- 
determined binary sources would not be significantly 
changed even if we do arrange follow-up observations tar- 
geted at the GW events. 

On the other hand, in the above estimate, it may be 
too simple to assume that the probability is indepen- 
dent of the mass and type of the galaxy and that GW 
events randomly occur in one of the galaxies. In fact, 
GW events are likely to occur in more massive luminous 
galaxies, and the fraction of redshift determination of bi- 
nary sources would be larger. In that sense, our estimate 
can be considered as the worst case and contains large 
uncertainty. 

Another way to determine the source redshift is to 
identify the electromagnetic counterpart when a NS bi- 
nary merges, which is believed to be the short gamma-ray 
burst (GRB). From GW observations in DECIGO band, 
we would be able to alert the electromagnetic telescopes 
to the event typically a few years before the merger. 
Given the half opening angle of a short GRB is « 10 deg, 
about 2% of the short GRBs point toward us. Even if we 
assume that only 10% of them have measurable redshift 
due to a noisy spectrum, dimming at high redshifts, and 
various telescope conditions, the fraction of redshift iden- 
tification is « 2 X 10"^, which is higher than that of the 
host galaxy identification. However, the opening angle 
of short GRBs has been poorly constrained [53| • Clarify- 
ing whether this method provides adequate counterparts 
needs further observation and more statistical samples. 

From the above estimations, it turns out that success- 
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ful rcdshift identification seems difficult and rather non- 
trivial. This indicates that the successful identification 
with high efficiency generally requires a dedicated follow- 
up mission with a sophisticated strategy. Therefore, the 
standard sirens without source redshifts guarantee mini- 
mally achievable FoM and strongly supports the feasibil- 
ity of the space-based GW detectors as a high-precision 
dark energy probe. 

Finally, note that the future galaxy surveys provide the 
redshifts of host galaxies below z « 2, while GW events 
distribute out to z « 5. However, from Eq. (U), most of 
GW events are localized around z = 1 — 2. In addition, 
the high-z GW sources are detected with less SNRs, and 
the existence of the dark energy would be less important 
at high-z. From these reasons, the follow-up observations 
at high-z, which is difficult and time-consuming, are not 
significantly crucial for the sensitivity of the standard 



where the overlap reduction function is set to unity j60| . 
The frequency range is set to /min = 0.2 Hz to escape 
from the WD confusion noise, and /max = 100 Hz as 
in the previous section. The power spectrum 5"^ is 
given in Eq. ([7]), but here it is multiplied by a factor 
(l/v^)^(-\/3/2)^ for the consistency of the definitions of 
the signal and noise in Eq. ((20| . The step function guar- 
antees figw > ^2gJ (/) 7?.NS , which is a necessary condi- 
tion for the clean subtraction of NS binaries. iVcorr is the 
number of independent correlation signals whose overlap 
reduction function is unity. Since DECIGO and BBO 
have two clusters at the same location and each cluster 
is composed of two independent (45°-rotated) interfer- 
ometers [6l[, iVcorr = 2. 



B. Results 



VI. PRIMORDIAL GRAVITATIONAL WAVES 

Next we will move on to the sensitivity of DECIGO 
and BBO to a primordial GW background. To search for 
the primordial GW background, one needs to accurately 
identify the waveform of NS binaries and subtract them 
up to the level sufficiently below the amplitude of the 
primordial GW background. In that sense, the parameter 
determination of the NS binaries and the detection of 
primordial GW background arc not independent topics 
but are closely related with each other. 



A. SNR formula 

Here wc consider a stochastic GW background that is 
(i) isotropic, (ii) stationary, (iii) Gaussian, and (iv) unpo- 
larizcd (see [5a] for the detailed discussions). To distin- 
guish the GW-background signal from stochastic detec- 
tor noise, one has to correlate signals between two detec- 
tors whose instrumental noises are uncorrelated with each 
other. The correlation-analysis technique has been well 
developed by several authors (ssl - fSTj and has been exten- 
sively studied in the detector configuration of DECIGO 
and BBO [l^ [5^, [s^. Under the assumptions that 
the energy density of the primordial GW background, 
r2gw(/), is independent of frequency and that the ampli- 
tude of its GW signal is smaller than that of instrumental 
noise, the SNR formula for two interferometers located 
at the same position, which is the case of DECIGO and 
BBO, is written as 



20\/2 7r2 



df 



1/2 



(20) 



Using Eq. ([^0]). we calculate ilg^r detectable with 
SNR = 5. Table IIVI and Fig. |S] summarize the resul- 
tant f2gw for various Tobs and ng . The sensitivity to J7gw 
in all cases with r„ = 1/2, 1/3, 1/5 are independent of 
the binary- merger rate no, because all NS binaries are 
subtracted out. As for the r„ = 1 case, some fraction of 
NS binaries, which is proportional to no, limits the sensi- 
tivity to rigw. For no = 10^^ Mpc~^ yi'^^, since a part of 
the noise curve is dominated by the NS-binary residual 
and the assumption that the amplitude of the GW back- 
ground is smaller than the noise is no longer valid, the 
sensitivity based on Eq. ([20]) is not shown in Table IIVI 
As explicitly understood from Eq. (|20l) , the sensitivity to 

i7gw is improved proportional to and T^^^J'^ except for 
the case with r„ = 1, where the sensitivity is limited by 
the NS-binary residual. 



VII. DISCUSSIONS 

A. Impact of GW observations on cosmology 

In the previous sections, we presented the sensitivity to 
dark energy and a GW background of a space-based de- 
tectors, which are summarized in Fig.|9l In what follows, 
we discuss some scientific consequences and the detector 
sensitivity. 

From cIl - X measurement, FoM for dark energy com- 
puted in Sec. IVl ranges from ~ 0.1 to ^ 10'^, depend- 
ing on the observation time, the merger rate, and the 
instrumental-noise level. Since DECIGO and BBO will 
be launched in the late 2020s, the FoM should be com- 
pared with other future projects of an electromagnetic 
observation probing dark energy at that time: type-la 
supernovae, baryon acoustic oscillation, or weak-lensing 
surveys. So typical criteria in FoM would be from 10 
to 100, which correspond to future projects of stage III 
and IV in the dark energy task force [5^, respectively. 
To achieve these criteria, observation time longer than 
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Tabs = 3 yr 


r„ 


ho = 10" 


V 




ho = 10" 


-e 


ho = 10" 


b 


1 


5.97 X 10" 


15 




5.97 X 10" 


-14 






1/2 


9.05 X 10" 


17 




9.05 X 10- 


-17 


9.05 X 10" 


-17 


1/3 


4.02 X 10" 


17 




4.02 X 10" 


-17 


4.02 X 10" 


-17 


1/5 


1.45 X 10" 


17 




1.45 X 10" 


-17 


1.45 X 10" 


-17 



Tohs = 5 yr 




ho = 10"' 


no = 10"" 


no = 10"'' 


1 

1/2 
1/3 
1/5 


5.94 X 10"''' 
7.01 X lO"''^ 
3.12 X lO"''^ 
1.12 X 10"" 


5.94 X 10"'* 
7.01 X 10"" 
3.12 X 10"" 
1.12 X 10"" 


7.01 X 10"" 
3.12 X 10"" 
1.12 X 10"" 




Tobs = 10 yr 


r„ 


ho = 10"' 


no = 10"" 


no = 10"" 


1 
1/2 
1/3 
1/5 


5.92 X 10""'' 
4.96 X 10"" 
2.20 X 10"" 

7.93 X 10"^* 


5.92 X 10"'* 
4.96 X 10"" 
2.20 X 10"" 

7.93 X 10"'** 


4.96 X 10"" 
2.20 X 10"" 
7.93 X 10"'* 



TABLE IV: Sensitivity to Slgw (SNR = 5). no is in the unit 
of Mpc"'^ yr"'. The symbol " — " in the table indicates that 
the NS residual partially covers the instrumental-noise curve 
and the SNR formula in Eq. (|20p cannot be used to compute 
the sensitivity. 
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FIG. 8: Jlgw detectable with SNR = 5 for each parameter 
set of r„, Tobs, and no- The observation time is denoted by 
colors: Tobs = 3yr (red, solid), 5yr (green, dotted), 10 yr 
(blue, dashed). For a fixed observation time, the merger rate 
only affects the r„ = 1 case: the upper and lower lines corre- 
spond to no = 10"®, and 10"^ Mpc"^yr"', respectively. For 
no = 10"® case, no data is plotted due to the same reason in 
Table HVl 



5 yr is preferable. Since the FoM rapidly improves being 
roughly proportional to ■, the observation time is a 
crucial factor. Given 10-yr observation and the typical 
rate of binary mergers, the FoM is ^ 100 and is compara- 
ble to a stage- IV project in the dark energy task force. It 
should be emphasized that this method requires no red- 



shift information of the GW sources and is completely in- 
dependent of any electromagnetic observation. Of course, 
if the GW observation data are combined with those of 
other electromagnetic observation, the FoM will be much 
improved as in the case of type-la supernovae. On the 
other hand, identifying the source redshifts strongly as- 
sists the cIl - X measurement. If the fraction of source- 
rcdshift determination a is larger than 10"* in the case 
with 5 yr-observation and no = 10"^ Mpc"'^ yr"^. the 
FoM is mainly determined by the binary sources with 
redshift information. The order-of-magnitude estimate 
of the fraction of redshift determination could be about 
10""'* in the worst case. Hence, we can conclude that the 
standard sirens without source redshifts guarantee mini- 
mally achievable FoM. 

As for a GW background, as shown in Table IIVI the 
space-based GW detectors such as DECIGO and BBO 
have enough sensitivity to directly probe a primordial 
GW background. On the other hand, cosmic microwave 
background (CMB)-polarization experiments j63 - l65| are 
also sensitive to the primordial GW background at the 
cosmological-horizon scale at present. Although the fre- 
quency range is far apart from that of the direct-detection 
experiment, the detector sensitivities should be compared 
with one another if one targets the detection of infla- 
tionary GWs. This is because the spectrum of GW- 
background energy density figw generated by single-field 
slow-roll inflation is predicted to be red-tilted j66l - l68t and 
is related to the tensor-to-scalar ratio defined at CMB 
scale via the evolution of the slow-roll parameters as 

f^gw(/*) = 4.81 X 10"i5V/r (^j^ . 

Here /q = 3.24 x 10"^* Hz, and is the frequency of the 
GW detector, e.g. 0.1 — IHz. tit is the spectral index 
of the tensor mode. Although the tilt and running of the 
spectrum have not been observationally well-constrained, 
the upper bound on the magnitude of the GW spectrum 
in the direct-detection frequency band is obtained in the 
limit of de-Sitter inflation (tit = 0) as 

^^gw(./*) = 4.81 X 10"^^ . 

Here we set /i72 = 1. Therefore, if the inflationary GWs 
are not detected by the CMB experiment, the GW de- 
tectors at least should have better sensitivity than the 
preceded GMB experiments for the first detection of 
an slow- roll inflationary GWs [7l|. Current constraint 
on the tensor-to-scalar is r < 0.2 [t^, which is trans- 
lated into figw < 9.6 X 10"is at / = /*. The future 
GMB-polarization experiments would be able to achieve 
r = 0.01, corresponding to Ogw = 4.8 x 10"^^ at / = 
To surpass the GMB-experiment sensitivities, r„ < 1/3 is 
required for the noise spectrum of the GW detector. Note 
that the current preconceptual design of the BBO noise 
curve marginally satisfies this sensitivity requirement. 




FIG. 9: FoM and figw achievable with GW observation alone (77 = 0) for each parameter set of r„, Tabs, and no. From the left 
to the right panels, the observation time is Syr, Syr, and 10 yr. In each panel, the merger rate is represented by curve types: 
no = 10~* (dashed), 10~^ (solid), lO"'^ (dotted), in the unit of Mpc~^yr~^. On each curve, points from the top to the bottom 
correspond to r„ = 1/2, 1/3, and 1/5. The horizontal lines correspond to r = 0.1 and 0.01 at CMB scale in the de-Sitter 
inflation. 



B. Other potential systematics 

In this paper so far, we consider systematic errors in 
the luminosity distance but not in the cosmological 
phase shift X. Since the cosmological phase shift is iden- 
tical to a cosmological redshift drift, the signal could be 
contaminated due to the peculiar acceleration of a source 
relative to an observer and the change of gravitational 
potential along the line of sight. According to Ref. [73j . 
such systematic noises are estimated, based on a linear 
perturbation theory. They found that the dominant con- 
tribution is the peculiar acceleration of the source, but 
the magnitude is smaller by one- or two- orders than the 
cosmological signal in the t ypi cal range of a redshift, i.e. 
z « 0.1 — 3. In Refs. [t^ l75l|. taking the nonlinear ef- 
fect of density perturbations into account, the authors 
considered a source in a galaxy and a galaxy cluster. In 
this case, the peculiar acceleration of the source is of the 
same order as the cosmological signal for a single source. 
However, the systematic error would be random for each 
source and can be reduced by averaging it out. Therefore, 
for our observation, we can conclude that the peculiar ac- 
celeration does not contribute to the cosmological phase 
shift of a GW. 

Another possible error that affects the sensitivity to 
a primordial GW background is subtraction noise that 
comes from parameter mismatches of each NS binary in 
the subtraction process. Although we neglect this con- 
tribution in this paper, it should be considered seriously 
because clean subtraction is crucial for detecting a pri- 
mordial GW background. Such a numerical investigation 
in the BBO detector configuration has been performed by 
Harms et al. [t^ . They showed that the residual noise is 
significantly suppressed by the subtraction process with 
the projection method and that BBO is able to detect the 
GW background with 0,^^ «2.5xl0~^''. Their result is 
obtained in the situation that ~ 10^ NS binaries are de- 
tected during 3-yr observation. Longer observation time 
would lead to the same conclusion, because the projected 



spectrum decrease with l/Tobs- However, they assumed 
the pessimistic merger rate, no = 10"'' Mpc~^ yr"-'^, 
which is 10 times smaller than the rate we used in this 
paper. If the rate is larger, the applicability of their re- 
sult should be reconsidered. Since the BBO sensitivity is 
comparable to DECIGO with r„ = 1/3, it is obscure how 
DECIGO with worse sensitivity affects the result. This 
issue of parameter mismatches should be investigated fur- 
ther in more broad parameter range in the future work. 



C. Signal overlaps 

DECIGO wiU detect ~ 10<^ of NS binaries, whose sig- 
nals could overlap and become individually indistinguish- 
able. To apply the projection method for the parameter- 
mismatch noise mentioned in the previous subsection, we 
at least need to subtract nearly matched signals. In this 
subsection, we discuss under what conditions these over- 
laps cause a problem in separately extracting each signal. 

The number of inspiral GW signals AN{f) in a bin of 
minimum frequency resolution A/ = 1/Tobs is given by 



A^(/) 



dN fdf 



dt \dt 



A/ 



(21) 



Here dN/dt is the merger rate per unit time, which is 
roughly 



no 



dN , 

dt V 10-6Mpc-^yr-i 



Substituting the frequency derivative [35|, |36[ 

^^^^8/3^5/3/11/3 
dt h ^ ■/ ' 

into Eq. (|2ip and setting AiV(/conf) = l, we obtain the 
critical frequency below which more than two signals are 
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in the same frequency bin as 



fci 



dN 



3/11 



^_8/ll^.-5/ll 



Roughly speaking, below this frequency, each signal is 
not distinguishable, which prevents us from extracting it 
individually. 

If we assume Tobs = 5 yr and that all signals are at z = 
1 for simplicity, the frequency /conf is 0.04 Hz for fiQ = 
10-'^Mpc"Vr"\ 0.07Hz for uq = 10"^ Mpc" 
and 0.14Hz for fiQ = 10^^ Mpc^"^ yr^^. As seen from 
Fig. [U the sensitivity to dark energy is hardly affected 
for the first two choices of no. However, if no = 
10~^ Mpc"'^ yr~^ and the lower cutoff frequency is set 
to /mill = 0.2 Hz, the FoMs are considerably deterio- 
rated from 9.2 to 6.6 x 10"'^ for r„ = 1 and from 81 to 
7.5 X 10~^ for r„ = 1/3. These results are expected from 

the scaling relation FoM oc T^^/*, since /i„in = 0.2 Hz 
corresponds to about 0.35 yr before the merger. To avoid 
this overlap problem, one needs more ingenious analysis 
method. We will leave it as a future work. As for the de- 
tection of an inflationary GW background, since we have 
already set /min = 0.2 Hz in our analysis, the signal over- 
lap does not directly affect the SNR. The degradation 
of the SNRs for NS binaries might increase the residual 
foreground though the parameter mismatches discussed 
in Sec. lVIIBl However, if they are projected out [t^, the 
residual would not affect the sensitivity. 



D. Comparison with other methods of the 
standard siren without redshifts 



There are two other proposals of the standard siren 
without electromagnetic counterparts, in which the de- 
generacy between a source redshift and a chirp mass is 
broken by observing a tidal effect on a GW |77l and by 
constraining the redshift range of each binary source with 
the observed chirp mass distribution of NSs [zll . In this 
subsection, we briefly comment on the comparison of our 
method with others. 

In the method utilizing the tidal effect [t^], we need 
to a priori know the equation of state of a NS. Further- 
more, the tidal correction on the GW phase is a 5 PN- 
order effect and can be observed only around the frequen- 
cies slightly before the merger, e.g. ^ kHz. Thus, this 
method would be available for Einstein Telescope (FT) 
[t^ and is completely complementary to our method. 
It will be interesting to investigate the multiplier effect 
when NS binaries are observed in both low and high fre- 
quency ranges. In the method with a chirp mass distri- 
bution [rsj , the sensitivity strongly depends on the width 
of the distribution. The authors found that in order to 
measure Hq with 10% accuracy by advanced LIGO 
[i^ . assuming ~ 100 sources detected, the half- width of 
the mass distribution should be within 0.04 Mq. This 
method is also applicable to FT and space-based detec- 
tors. It is valuable to compare the sensitivities of these 



methods with that of our method, but it is beyond the 
scope of this paper and we leave it as a future work. 



VIII. CONCLUSION 

Proposed space-based GW detectors, DECIGO and 
BBO, can be a unique probe for dark energy and a pri- 
mordial GW background. We investigated these target 
sciences from integrated point of view, allowing their 
noise curves to be scaled appropriately and including 
confusion noises produced by astrophysical sources. Us- 
ing millions of NS binaries detected during the obser- 
vation, we estimate the sensitivity to the parameters of 
the equation of state of the dark energy with/without 
identifying the redshifts of host galaxies. As a result, 
we found that the detectors without the redshift infor- 
mation of the sources have a constraining power com- 
petitive to the future electromagnetic observations. This 
is a great advantage since the GW detector alone can 
probe for the cosmological expansion and enables us to 
compare the data of purely GW observation with those 
obtained in other electromagnetic observations. With the 
help of the redshift information, FoM « 100 correspond- 
ing to the stage HI or IV of the dark energy task force is 
easily achieved with at most a few hundreds of sources. 
The standard sirens without redshifts guarantee mini- 
mally achievable FoM and strongly supports the feasibil- 
ity of the space-based GW detectors as a high-precision 
dark energy probe. The detection and cleaning of GWs 
from the NS binaries are also necessary to search for a 
GW background generated in the early universe. We 
took the subtraction procedure of the NS binaries into 
account and computed the sensitivity to the amplitude 
of a GW background. As a result, the default BBO sensi- 
tivity or the DECIGO with three-times-better sensitivity 
is marginal in order to successfully subtract the NS fore- 
ground and to surpass the future CMB-experiment sensi- 
tivity to a primordial GW background. In summary, the 
subtraction of NS foreground is significantly important 
for both dark energy and primordial GW background, 
and our study here provides useful information to opti- 
mize the experiment. 
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Appendix A: Fitting formula for NS residual 

The NS residual fraction 7?. ns as a function of pth/ Pc(5) 
is given by Yagi and Scto in [29| . For the convenience of 



the analysis in this paper, we fit their result and provide 
the fitting formula |8l|. Fitting TZ^s by polynomials up 
to 5th order in x = pth/Pe(5) leads to the following piece- 
wise expression: 



f -3.695 + 2.297 x 10^ x ~ 5.714 x 10^ + 7.105 x 10^ x^ - AAU x 10^ x"^ + 1.098 x 10^ x^ 

for X < 0.83 

9.675 X 10-2 - 5.663 x IQ-i x + 1.321 x^ - 1.533 x^ + 8.847 x lO^^ x^ - 2.029 x IQ-i x^ 

for 0.83 < a; < 0.9 



2.166 X 10^2 _ 1^ 255 x IQ-i x + 2.815 x lO^^ x 



NS 



1 -2 _ 3.039 X IQ-i x^ + 1.563 x 10"! x^ ~ 2.993 x lO'^ 

for 0.9 < x < 1.1 



3.561 X lO-'' - 5.516 x 10-^ x + 2.321 x 10"^ x^ - 3.741 x 10"^ x^ + 2.279 x 10"^ x'^ - 4.424 x lO'^ x^ 

for 1.1 < x < 1.8 

2.594 X 10-2 - 8.338 x lO'^ x + 5.340 x lO'^ x^ - 6.790 x 10"^ x^ + 3.560 x lO""* x^ + 6.722 x 10"^ x^ 

for 1.8 < X 



Appendix B: Comparing between error sizes of 
and X 



In this Appendix, we estimate the error size of the lu- 
minosity distance in the presence of possible system- 
atic errors and then compare it with that of cosmological 
phase correction X. 

The error size of the luminosity distance can be esti- 
mated via the same procedure as in Sec. IV B] but includ- 
ing systematic errors. It is known that the luminosity 
of a GW is magnified by the matter inhomogeneities of 
large-scale structure along the line of sight (e.g., [H-d^), 
which systematically changes the luminosity distance to 
each binary system and contributes to the measured lu- 
minosity distance as a systematic error [s^. In addition, 
the peculiar velocity of the binary along the line of sight 
randomly contributes to measurement error via Dopplcr 
effect [87| . These systematic errors to the luminosity dis- 
tance are summarized as 



AdLiz) 



dL{z) 



t^inst(2)+Criens(2)+'^pv(^) : 



with 



0-\ens{z) = 0.066 



l-{l + z) 



-0.25 



n 1.8 



O-pv(z) 



1 



0.25 



(Bl) 



(B2) 



H{z)dL{z) 



v.gal 5 



where (Xinst, ciens, and CTpv are induced by the instru- 
mental noise, the lensing magnification, and the pecu- 
liar velocity of binaries, respectively. Cv.gai is the one- 
dimensional velocity dispersion of the galaxy. Taking into 



■2 10-^ 

T3 



Ole 




instrumental, Syr 
instrumental, Syr 
■ instrumental, lOyr 
lensing 
peculiar velocity 



Op 



z 



FIG. 10: Measurement accuracy of the luminosity distance 
with a single binary as a function of redshifts. The curves 
tagged (Jinst are those determined only by instrumental noise 
and with the observation time Syr (red, solid curve), Syr 
(green, dotted curve), and 10 yr (blue, short-dashed curve), 
respectively. The lensing error and the peculiar velocity error 
are represented by magenta (long-dashed) and light blue (dot- 
dashed) curves. 



account the nonlinear effect of gravity, it is often set to 
aVgai = 300km s~^, mostly independent of the redshifts 
[88| . In Fig. [TUl the errors of the luminosity distance 
are illustrated. The lensing error dominates at almost all 
redshift range. 

We present Aiuq, Awq, and FoM in the case with a = 1 
in Table [V] As anticipated, z - measurement is much 
more powerful than the dL - X measurement. 

Next, we show that the error of d^ is much smaller than 
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Tohs = 3 yr 




no = 10"' 


no = 10"" 


no = 10"" 


r„ 




AWa 


FoM 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


1 
1/2 

1/3 
1/5 


7.833 X 10"* 
5.202 X lO"'' 
4.537 X 10"* 
4.150 X lO"* 


1.254 X 10"^ 
8.782 X 10"^ 
7.828 X 10"^ 
7.265 X 10"^ 


1.101 X 10" 
2.272 X 10^ 
2.877 X 10^ 
3.357 X 10^ 


2.506 X 10"* 
1.692 X 10"* 
1.435 X 10"* 
1.312 X 10"* 


4.007 X 10"^ 
2.845 X 10"=^ 
2.475 X 10"=* 
2.297 X 10"^ 


1.078 X 10" 
2.162 X 10" 
2.877 X 10" 
3.358 X 10" 


8.696 X 10"" 
8.782 X 10"* 
4.537 X 10"" 
4.150 X 10"" 


1.378 X lO"'' 
4.632 X lO"'^ 
7.625 X 10"* 
7.265 X 10"* 


9.095 X 10" 
2.272 X lO'^ 
2.877 X lO'^ 
3.358 X 10'^ 




Tabs = 5 yr 




no = 10"' 


no = 10"" 


no = 10"" 


r„ 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


1 
1/2 

1/3 
1/5 


5.980 X 10"* 
3.997 X 10"" 
3.498 X 10"* 
3.213 X 10"* 


9.611 X 10"-' 
6.762 X 10"^ 
6.042 X 10"^ 
5.619 X 10"^ 


1.875 X 10'' 
3.833 X 10^ 
4.829 X 10^ 
5.613 X 10^ 


1.913 X 10"* 
1.264 X 10"* 
1.106 X 10"* 
1.015 X 10~* 


3.071 X 10"^ 
2.138 X 10"^ 
1.911 X 10"^ 
1.777 X 10"^ 


1.836 X 10" 
3.833 X 10" 
4.829 X lO" 
5.613 X 10" 


6.638 X 10"" 
3.997 X 10"" 
3.498 X 10"" 
3.208 X 10"" 


1.056 X 10"-' 
6.762 X 10"* 
6.042 X 10"* 
5.619 X 10"* 


1.549 X 10' 
3.833 X lO'^ 
4.829 X lO'^ 
5.613 X lO'^ 




Tobs = 10 yr 




no = 10"' 


no = 10"" 


no = 10"" 


Tn 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


Awo 


AWa 


FoM 


1 

1/2 
1/3 
1/5 


4.187 X 10"* 
2.813 X 10"* 
2.468 X 10~* 
2.267 X 10"* 


6.746 X 10"-' 
4.764 X 10"^ 
4.265 X 10"^ 
3.971 X 10"^ 


3.803 X lO'' 
7.718 X 10^ 
9.691 X 10^ 
1.124 X lO'^ 


1.339 X 10"* 
8.895 X 10"^ 
7.805 X 10"^ 
7.170 X 10"^ 


2.155 X lO""* 
1.507 X 10"^ 
1.349 X 10"^ 
1.256 X 10"^ 


3.725 X 10" 
7.718 X 10" 
9.691 X 10" 
1.124 X 10" 


4.642 X 10"" 
2.813 X 10"" 
2.468 X 10"" 
2.267 X 10"" 


7.405 X 10"* 
4.764 X 10"* 
4.265 X 10"* 
3.971 X 10"* 


3.148 X 10' 
7.718 X lO'^ 
9.691 X lO'^ 
1.124 X lO'^ 



TABLE V: Measurement accuracy of equation of state of dark energy and FoM when assumed that redshifts of all binaries are 
identified by electromagnetic follow-up observation of host galaxies {a = 1). no is in the unit of Mpc"^ yr"^- 




FIG. 11: Errors of X{z)/Ho as a function of a redshift. The 
Solid curve denote the error due to instrumental noise in a 
GW observation, and the dotted curve denote the error prop- 
agating from the lensing error of d_L. Assumed parameters are 
T-n = 1, Tobs = 3yr, and no = 10"" Mpc"^ yr"\ 



tion of a redshift z. DifFcrcntiating Eqs. and 
with respect to z and eliminating z give 



AX{z) 



where 

dH{z) 
dz 



H-{\ + z)^Ml 



2(1 + z) i+(i±£)! 



1 



(B3) 



{l + z)H{z) 



l+WQ+Wa 



1 



Q{z) + -H'^^m{l + zY 



Q[z) = -Hlil - a„)(l + 



X cxp 



l + z 



that of X and that can be regarded as a corresponding 
redshift in the fiducial model of the universe when we 
constrain cosmological parameters with standard sirens. 

To compare the error sizes of dh and X, let us see 
whether the error of affects uncertainty of X or 
not. For simplicity, we fix cosmological parameters to 
the fiducial values: wq = —1, Wa = 0, il„i = 0.3, 
i?o = 72kms"^ Mpc"^ and regard d^ and X as a fimc- 



Dominant error contribution to d^ is the lensing error, 
explicitly given in Eq. (|B2[) . After substituting uions 
for Ad^/dL, we plot numerically evaluating AX{z) in 
Fig. 1111 assuming parameters = 1, Tobs = 3yr, and 
TT-o = 10"^ Mpc"^ yr"^, and compare with AX{z) com- 
ing from instrumental noise in a GW observation. We 
found that the error of is much smaller than that of 
X and can be neglected. This conclusion does not change 
if other parameters for r„, Tobs, and no a-rc used. 
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